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COMPUTATION OP SUBSONIC FLOV/ AROUND AIRFOIL SYSTEMS 
WITH MULTIPLE SEPARATION* 

Klaus Jacob** 

SUMMARY 

A numerical method for computing the subsonic flow around 
multi-element airfoil systems has been developed, allowing for 
flow separation at one or more elements. Besides multiple rear 
separation also short bubbles on the upper surface and cove bubbles 
'v can approximately be taken into account. Also, compressibility 
effects for pure subsonic flow are approximately accounted for. 

After presentation the method is applied to several examples and 
improved in some details. Finally, the present limitations and 
desirable extensions are discussed. 

NOTATION . ***/£ 

superposition factors for the corresponding basic solutions 
) > C3-2 ) 

m-th connection point for the dead water on profile 
number 1 

superposition factor for basic solution (b^), i.e., pure 
circulation flow around profile number 1 
lift coefficient 

pressure coefficient ■ (p- p^)/(v*p/ 2 ) 

pitch moment coefficient, referred to the 1/4 point 
superposition factor for the basic solution (d^ ^) , i.e., 
pure outflow from profile number 1, where the source occu- 
pation starts at point 
yield of source occupation 
trailing edge 

* DPVLR Research Report 81-24, 1981, 63 pages. 

Institute for Theoretical Fluid Mechanics of DFVLR, Goettingen. 

*** 

Numbers in margin indicate pagination of foreign text. 
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geometric functions defined by equation (2d) and equation 
(2e) 

profile chord 

Mach number of Incident flow 

number of given support point on the contour of profile 
number k ‘ • 

number of profiles In the profile system 

point on the upper boundary line of the dead water region 
of profile 1 

static pressure and static pressure of Incident flow 
respectively 

Reynolds number = V^*l kinematic viscosity 
running length along profile contour or gap distance 
(In Figure 13) 

various lengths In Figures 4 and 13 

attachment point of dead water region along the underside 
of profile 1 
local flow speed 
Incident flow speed 

Incident flow speed of resulting flow 

normal and tangential component of flow speed along profile 
contour. These are Identical with the intensity of a sur- 
face source and vortex occupation 
X and y components of flow speed 
cartesian coordinates 
angle of attack 

displacement thickness and momentum loss thickness of 
boundary layer 

parameters; to each value of for example there Is the 
unique contour point of the profile number k associated 
with it 

circulation around profile 1 
adiabatic exponent 

transformed source Intensity according to equation (2b) 
transformed vortex Intensity according to equation (2a) 
density of flowing medium 
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Subscripts : 

i, j, k running numbers of profile 

m, y, V running numbers of points on profile contours 
<» in the flow far away from the system 

LA at the separation point of the laminar boundary layer 

1. INTRODUCTION 

Wings with slats and slot flaps (mechanical lift aids) have 
for a long time been used te improve the takeoff and landing per- 
^ formances of aircraft. The longitudinal section through such a 

wing system results in a profile system (See, for example Figure l4, 
top). 

Over the last few decades intensive experimental and theoret- 
ical profile research has led to substantial improvement of indi- 
vidual profiles, for example, substantially reduced drag for fast 
flight due to the development of laminar profiles and over-critical 
profiles. There has also been improved high lift and stall behavior 
by suitable curvature, S-shape, etc. Improvements to the (two- 
dimensional) profiles have also resulted in better (three-dimen- 
sional) wings especially for large wing aspect ratios. The calcula- 
tion possibilities for individual profiles have advanced consider- 
ably. For incompressible, frictionless flow, very exact surface 
singularity methods are available. Friction and even flow separa- 
tion can be approximately considered in calculations- 

On the other hand, for profile systems, the magnitude of avail- 
able experimental data as well as the possibility of calculation 
/• are rather restricted. The high lift configurations which are very 
Important for takeoff and landing of our aircraft have developed 
empirically for the most part. The calculation of such profile 
systems can be done for frictionless flow using singularity methods 
(see, for example, [1-3]), but friction and partial separation are 
even more important in practice than is the case for individual 


3 



profiles. Calculation methods for profile systems with friction 
and separation are still in the beginning stages ’of development. 

The aircraft industry is very interested in reliable and practical 
calculation methods of this type in order to develop better high 
lift configurations for new aircraft and to reduce the experimental 
expenditure. 

Based on the calculation method for the individual profile with 
flow separation developed by the author, during the 70 ’s 
several similar methods were developed for profile systems, for 
example, z. B. von BHATELEY and BRADLEY [5], JACOB and STEINBACH 
[6], HAYASHI and ENDO [?] and OLSON [8]. All of these methods have 
one common feature in that a potential flow with an imbedded simu- 
lated dead water region is calculated using surface singularities 
and this is combined with boundary layer calculations using itera- 
tions. Differences between the individual methods occur in the 
modeling of the dead water region and in the boundary layer method 
used. All of the methods up to the present were restricted to anal- 
yzing pure trailing edge separation^ and additional possible physi- 
cal phenomena such as separation bubbles, bottom side separation 
compressibility effect were not considered. Wake boundary layer 
confluence is only included approximately in [8]. A ground influ- 
ence is only considered in [6,6b]. Most methods are restricted to 
simple separation, that is, separation at only one element of the 
profile system. 

A general and really reliable calculation method can only be 
developed in steps by researching all of the possible physical phen- 
omena, continuously comparing them with experiments, theoretical 
modeling and these all have to be included in the calculations. In 
addition, reliable criteria for the occurrence of individual pheno- 
mena are required. Of course, a restricted method can be useful 
but only if its limitations are known and it is appropriately 
applied . 





In the present paper we will extend our method for profile 
systems [6] by three steps: 

The most important part of the work was the extension to cases 
of multiple trailing edge separation. This means that simultaneous- 
ly separation regions are allowed at several parts of the profile 
system, and their interaction is considered. It is only in this 
way that an exact determination of the maximum lift, and calcula- 
tion of the lift curve in the post-stall range becomes possible. 

In addition, two smaller extensions are made which allows 
the consideration of short separation bubbles and a consideration 
of the Influence of compressibility for subsonic flow. /II 

The extended method is then tested with various examples and 
new application limits and desirable further extensions are deter- 
mined. 

■j 

2. THE CALCULATION METHOD 


2 . 1 Introduction into the previous method [6,6a] 

We will first consider the simplest case of a flow around a 
profile system with separation, that is, incompressible flow with 
a simple trailing edge separation with a large Reynolds number (Fig- 
ure 1). It is assumed thcit only three of the basic elements of the 
flow field exist, that is, approximately frictionless external 
flow, the boundary layer at the body contour and a dead water 
region with almost quiet flow, above add behind one of the profiles. 
First of all, no additional elements are assumed (for example, 
separation bubbles and further dead water regions) which could sub- 
stantially change the flow. 

The calculation of such a flow requires the calculation or 
theoretical modeling of each of the elem.ents and consideration of 
their interactions. 
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Our calculation Is divided into two main parts which are 
coupled by Iteration: 

1) the calculation of a potential flow with enclosed simulated 
dead water displacement region for a given separation point, 

2) boundary layer calculation at the body contour outside of the 
dead water region* for a given pressure distribution. 

The coupling is carried out through the pressure distribution 
and the separation point: 

A separation point is found by iteration for which the corres- 
ponding potential flow has a pressure distribution where the cor- 
responding boundary layer separates exactly at this point. 

Calculation of the potential flow with dead water : 

For the calculation of the potential flow, a singularity method 
with vortex occupation on the body surface is used. In addition, 
on the profile with separation a source occupation is applied through 
which an outflow region is generated which is similar to a dead 
water region in its characteristics regarding the displacement 
effect with respect to the external flow. For this simulated dead 
water region, it is required that it has constant pressure along 
its boundaries near the profile. 

After specification of the source occupation v^(<f>), there is 
first the problem of determining the vortex occupation in such 

a way that the kinematic flow condition at the body contour is sat- 
isfied everywhere. For solving this problem, for the case of a 
single profile of MARTENSEN and SENGBUSCH [9], the following inte- 
gral equation of the second kind was derived: 



For the case of a profile system consisting of N profiles, equa- 
tion (1) can be extended to the;)following system of coupled linear 
integral equations 


( 2 ) 




N 2n 


j = 1 0 




2it 

= 2v^(x(»j^)cosa +y(ipj^)sin a) + L ( d* 

0 


for k = 1 ,2, . . .,N . 


Here we have: 

i number of the profile with separation, along which there is the 
transformed source occupation w , 


j,k running profile numbers 


(2a) 


(2b) 

•x)n('ii) = 

(2c) 

s (®) = ds/d(p = 

= /(dx/d®)»+ (dy/d<p)2 


transformed vortex intensity at 
the point (}ij^ of the profile k 

transformed source Intensity at 
the point of the profile i 

derivative of the running length 
s along the profile contour with 
respect to the parameter <j) 


It is assumed that each profile is available in a parameter 
representation x(«), y(ip) so that starting with'” = o at the trailing 
edge, the profile is passed through exactly once in the counter- 
clockwise direction for "> = o to J • 

In equation (2), in addition and f are the Influence 

functions. They are the tangential speed which is produced by unit 
vortex or unit source which lie at a point defined by ip, at a con- 
tour point defined by 4>. It is calculated according to the following 
formulas: 
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(2d) 

(2e) 


y(C)^) • Ix(ipj^) - ] - x(oj^) [y(<i)^) - 1 

[x(9j^)-x(i^j)]’ + [y(9j^)-y('l’j)l' 

x(9^) : [x(«j^) - x(i^) 1 + y(»^) • ty (®^) - yU^) 1 

[x(®j^)-x(tjij^)]^ + [y(»j^) -y(i'j^)]’ 


» 

The system of Integral equation (2) can be transformed by 
discretization into a set of N linear equation systems 
which can be solved numerically. One condition for a unique 
solution is that not only all of the geometric data but also 
the incident data and a as well as the source occupation on 
the right side of equation (2) are known. . In addition, one addi- 
tional condition is required for specifying the circulation for 
each profile. 

By the simultaneous solution of the equation system for differ- 
ent right sides and different additional conditions, the so-called 
basic solutions given in Figure 2 are calculated and with a weighted 
superposition we finally are able to generate flows of the kind 
shown in Figure 1. /l^ 

The basic solutions (a^) and (a 2 ) (parallel basic flows) are 
obtained for v^ = 1, a = 0° and 90*^, = 0 respectively and the 

additional conditions ° for all k. ( ° means the 

trailing edge of profile k). Basic solution (b) (pure circulation 
flow around profile 1) is obtained from 'v„= o, o addi- 

tional conditions 

2n 

’’i ° / ® and =0 for all k 1. 

0 

One basic solution (d) (pure outflow) is obtained from v^ = 0, 
a source occupation different from 0 * on profile i ° 

between A and U and the additional conditions 

2n 

’’i = / " ° and oj^(»k=o)=o for- all k 1. 

0 

in ~ 

See [10] in section 2.2.1 for the selection of the source occupation. 
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The basic solutions of the last kind are made available for differ- 
ent positions of A = A , m = 1.2,... so that it is possible to pro- 
duce displacement regions of various length by superposition with 
other basic solutions. The basic solution (d) corresponding to 

A will be called (d ). 
m p 

After obtaining the basic solutions it is then necessary to 
produce a resulting flow by superposition with suitable, weighted 
factors a^, b and d^^^ which have a given incident flow V^, a 

and a displacement region along whose boundary the pressure is 
constant. 

The incident flow is only determined by the basic solutions 
(a^^) and (a2) which are the only ones not equal to zero. For 
= 1 (reference speed for all speed magnitudes) we find the following 
for the superposition factors a^ and a 2 ^ 

( 3 ) ^ ^2 ~ 

The remaining superposition factors b and d of the circulation 
flow and outflow essentially determine the direction and thickness 
of the displacement region. They are determined by iteration in 
such a manner that the resulting flow has equal pressure at three 
characteristic points, that is, at the two separation points A^ and /1 5 
U and a point 0 on the upper boundary of the displacement region 
above the profile trailing edge (see Figure 1). In this way we have 
satisfied the condition for constant pressure along the boundaries 
of the displacement region near the profile approximately. Therefore, 
one has to solve the two equations 

( 4 ) 

We will discuss this problem further in section 2.2.2. 
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Calculation of the boundary layer 


After we have obtained a potential flow with a simulated dead 

water region with a specified attachment point A^, we must examine 

whether this flow Is physically realistic. It Is only when the 

flow has a boundary layer which separates exactly at A^ that dt can 

really exist. Otherwise, the calculation for different A has to 

m 

be repeated until the correct solution Is found. During this Iter- 
ation process we require boundary layer calculations which will pro- 
vide reliable solutions which are not too complex. We use the Inte- 
gral method of Rotta In [11] which was programmed. 

Our* boundary layer calculation always starts lamlnarly at the 
front stagnation point and can be continued In a turbulent manner 
as desired starting at a specified transition point or according to 
the empirical transition criterion of Michael [10] or It can start 
from a calculated laminar separation point and this Is continued up 
to turbulent separation. In addition to determining the separation 
point, the boundary layer calculation Is also used to determine the 
friction drag and the boundary layer displacement thickness 6^. The 
boundary layer displacement In our calculation Is considered by add- 
ing an additional basic solution C^q) with the source condition 

(5a) (5a) = d(6j • v^)/ds, 

/1 6 

and the additional condition 

(5b) 

* 


'Recently, we have preferred the additional condition ■ a*) 

with a<is2i/3o, that Is, the condition of equal velocity at two 
points which lie near the trailing edge at the same distances 
along the top side of the bottom side. 
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In later sections we will discuss some new features In fur- 
ther details of the boundary layer calculation. ’ 

The overall procedure for the previous method 

The overall probedure includes a preliminary phase consisting 
of two Iterations inside of one another and a final calculation: 

The preliminary calculation contains: 

The preparation of all geometric quantities, calculation of 
basic solutions a^^ and a .23 

calculation of the attached frictionless flow (superposition of 
basic solutions (a^) and (a 2 ) and the corresponding boundary layer 
and finally the calculation of the basic solutions (b), (d^), (d^), 
(d2)... 

The internal iteration includes: 

The calculation of a flow with simulated dead water region 
for an attachment point at the profile number 1. This is done 
by superposition of basic solutions -(ai), (aj), (b), (do) and ' 

The determination of suitable superposition factors b and d^^ is 
done by iteration. The purpose of this iteration is to achieve 
equal pressure at points A , U and 0. 

I 

The external iteration is used to find a flow with simulated 
dead water region where the attachment point A^^^ coincides with the 
separation point A of the boundary layer. In addition to the inner 
iteration it includes a boundary layer calculation for the top side ' 
of profile 1 for determining the separation point A. Starting with 
an attachment point near the trailing edge, this is done in small 
steps and pushed forward until the boundary layer calculation for 
the first time results in no separation ahead of A^. 

The final calculation contains the following: 
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Dfe termination of the force and moment coefficients of the 
profile system by integration of the pressure and shear stress 
distribution, also boundary layer calculations for the top side 
and the bottom side of all of the profiles participating in the 
system. We then find whether as assumed separation actually /17 

occurs for a single profile (number 1) for the given configuration. 
Incident flow angle and Reynolds nuftiber. Otherwise, the calculated 
flow cannot be looked upon as realistic. 

2 . 2 Determination of multiple trailing edge separation 

For maximum lift and especially in the post stall range, pro- 
file systems often have several separations in regions at various 
parts of the system. For example, a profile with a highly deflected 
gap flap will only exceed its maximum lift at an angle of attack 
for which there is a substantial separation over the main wing for 
flap separation. In order to allow one to calculate the maximum 
lift and the post stall behavior for profile systems, it is necess- 
ary to expand our method to the case of multiple separation. 

In order to be able to not only simulate a profile but in the 
extreme case to simulate all profiles of the system and their 
separation regions, suitable source conditions along all profiles 
are required. This leads to the following expanded Integral equa- 
tion system 


N 2< 
j«1 0 

for '^k « 1,2,...N 

•where 

(6a) “ 2v^(x(»j^) cos a + y(*j^)slno) * 

N 2» 

* -hT, I 

j=1 « 
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As will be discussed in the following sections, this system 
of coupled linear integral equations allows one to calculate an 
additional supply of basic solution's which then can be superimposed 
during an additional iteration method in such a manner that the 
resulting flow contains several simulated dead water regions. /1 8 

2.2.1 The required basic solutions and their calculation 

In addition to the basic solution (a-j^) and {B. 2 ) discussed in 
section 2.1, for each profile i (i = 1.2,...N) we require one cir- 
culation flow (bj^) and a set of outflows (d^ ^), (d^^ 2 ) 9 * 

various attachment points A. ,, A. -,...of the source occupation. 

1 , ± 1 3 ^ 

We will now exactly define these basic solutions and give the 
corresponding right sides *‘*'*’k* for equation (6) as well as the 
additional conditions: 

Basic solution (b^): Source free, circulation flow which 

vanishes at infinity with specified circulation = 2 n 1 v„ around 
profile i, which satisfies the Kutta outflow condition at the trail- 
ing edge at all other profiles. This means: 
u)n(^pj). = 0 for all j 
= 0 

“211 . 

0 

“t = 0) = 0 for all k 1. 

The first and second conditions substituted in equation (6a) 

give 

R(<J)j^) = 0. 

The third and fourth conditions are additional conditions for unique 
calculation of these basic solutions. For each profile, a basic 
solution of this kind is calculated. 

Basic solution (d^^ Outflow which vanishes at infinity which 
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is caused by. a specified source condition between contour points 

* 

A. and U. has no circulation around profile i and satisfies the /19 

X j m X 

Kutta outflow condition at all other profiles. This means: 

Vco = 0 

u (.(. ) ^0 between A. and U. , w (if;.) = 0 for j i 

2it 

J = 0 

0 

u (»^ = 0 ) =0 for all k / i. 

t k 

The first and second conditions in equation (6a) give 

2tt 

0 

The source occupation which occurs here is calculated accord- 

ing to equation (2b) and the normal velocity distribution v^ accord- 
ing to Figure 3 is selected (see also [10]). 


The third and fourth conditions are again the additional con- 
ditions for unique calculation of this basic solution. For each 
profile, several basic solutions of this kind are calculated (for 
several attachment points of ^) on the profile top side. 

Numerical solution of the Integral equation system 


The integral equation ' system (6) after discretization* gives a 
linear equation system for each value of k. If is the number of 
support points ^k,u' “ ° qj^ profile k, then the k-th 

equation system contains equations. All of the N equation systems 

then give 


N 

” “ 2 *^k equations 

k=1 


Approximation of the Integral by sums according to the rectangle 
rule . 


in 



at all 


The number of unknown values of vortex Intensity 
support points ^ ° j = 1 , 2 ,. ..,n j^g 

N 

M = S “j • 

j = 1 

Nevertheless, it is not immediately possible to obtain any gefneral 
unique solution of the overall system because each Individual equa- 
tion system is simply linearly dependent. It is only after one addi- 
tional condition is added to each of the equation systems that one 
obtains a unique overall solution. (The summing up of the additional 
equations for all of the equations of the system has proven itself). 

In order to obtain good numerical results with the minimum 
effort, in earlier work a number of special measures for establish- 
ing equation systems were developed, for example, a method for 
obtaining favorable support point distribution [12], special formu- 
las for better treatment of relatively thin profiles and profiles 
with corner trailing edges [13] and many other things. In [10], 
there is a discussion of this. We will not discuss these details 
again even though they are important for practical success of the 
method. 

By solving the overall system of linear equations for those 
right sides and those additional conditions which correspond to a 
special basic solution (g), one finally obtains all of the values 
of the transformed vortex occupation • and, therefore, using 

equation (2a) the tangential speeds ^.ll of 

the support points ^jv' '' ° of all profiles j = 1.2,...,N for 

this special solution (g). 

In the practical calculation of the tangential speed for our 
different basic solutions, it is economical to simultaneously solve 
the linear equation system for all basic solutions in which the 
equation systems only differ by their right sides (including the 
additional conditions). This applies for the basic solutions (a^) 
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and (a 2 ) as well as for e^ch group (b^), (d^ ^), (d^^ 
a fixed i. 


2.2.2 Superposition of basic solutions 


/21 


After we have made available vortex and velocity distribu- 
tion for all of our basic flows. It Is now necessary by means of a 
suitable weighted supersposltlon to obtain a resulting flow with 
simulated separation regions at several profiles, and an equal pressure 
condition has to be satisfied for each of the separation regions. 


The result of a superposition 


( 7 ) 


(a ) (aj) 


(d„) 


N 


i=1 




depends on the values of the basic solutions and on the selected 
superposition factors ® 2 ' ^ 2 ' • 


For = 1 and a given outflow angle a, a^ and a 2 are to be 
looked upon again as given according to equation (3). 


For given attachment points of the separation regions ^ , 

1 = 1.2,...,N the 2 N additional superposition factor and ^ 

d. ,1=1.2,...,N are to be found In such a manner that the 
l,m.’ 

pressure coefficients 


( 8 ) 


Cp= 1-(v/v_^)^ - f(bj .dj ,bj ,dj 


are the same at the comparison points ^ and 0^ of each sep- 

aration region. In other words, the commoA zero of the two N 
functions 


( 9 ) 


'pA 







i = 1,2 


,N 
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is to be determined. 


The values of the function can be easily 'calculated for a 
given set of superposition factors using equations (7) and (8). 

The calculation of the function values of h^ is more complicated. 

First of all, the field points 0^ along the upper boundaries of the 
separation regions have to be determined and the velocities there 
v« have to be calculated. /22 

The velocity components v and v at an arbitrary field point 

X y 

P(x,y) are calculated from the following formulas: 


(10a) 


M. 


v^(x,y) - V^’cosa 




(lOb) 


v^U.y, = V.. sina . ^ ^ ^ . 


A field point 0^ can now be found by starting at a point near to 

A. and following the speed direction (streamline), until one has 
1 9 m . 

reached the desired position above the profile trailing edge. Then 

2 a a 

one obtains Vq from ^ ^y • 

Now that we are in the position of numerically calculating the 
values of the functions g^^ and h^^^ according to equation (9) for a 
given set of superposition factors b, ,...,d.. , we will now consl- 

der the determination of the zeros. 


We have to find a set of superposition factors so that 


(11) " ‘ 1 = 1,2 N 


h. (b, ,d , , . . /d., ) - 0 


' l,m 


N' N,mi 


‘N 


The solution of this problem is done using an iteration and recourse 
to physical interpretation. First of all^ we assume that 
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reasonable zero solutions are available for the superposition fac- 
tors bl"*, d j,. . . , and we will discuss this at the end 

of this section. 


We start the entire Iteration calculation with small separation 
regions at all profiles, that Is, with attachment points ^ which 
are only a few per cent of the profile chords away from the trailing 
edges . In this case, one can expect that the mutual Influencing 
of the separation regions Is small and, therefore, the correct value 
pair b^, d^ ^ will barely depend on the other ones. /23 


Each value pair b^, d^ ^ Is then calculated using a two-dimen- 
sional Newton method (see [ 6 a], and for all of the other superposi- 
tion factors, the previously highest approximations are used and 
fixed: o 


I and the highest approximations of tij,dj ^ for j 7 ^ 1 


-»b 


(n+i) ^(n+i) 


n = 0,1 »2, . 


As soon as this Iteration Is finished, the next profile Is cal- 
culated: 1 = 1.2,...,N. The result of this Iteration are Improved 
approximations ^ 1 ' ‘^ 1 , 1 ' • • • ' ‘^,1 of all of the superposition fac- 
tors. But Ih the calculation of, for example, ^ 1 ' » ^ only the 

zero approximation of the other superposition factors were available. 
"A second run through this Iteration may be useful (after substitu- 
tion of the zero approximations with new values) (especially If the 
separation regions are large and the zero approximations were not 
very good) . 

Using the selected superposition factors, the resulting over- 
all flow Is known. If It Is necessary to make calculations for 

additional attachment points A. (that Is. larger separation regions) 

1 , m . X 

(see section 2 . 2 . 3 ) j then this can be done In the same way as for 

We selected the first attachment point A. ^ as well as the attach- 
ment point U, on the lower side at about^’"^ 98.55? of the profile 
chord . 
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the attachment points One requirement Is that we will always 

have a good zero approximation for superposition ' factors 



This Is done by allowing the attachment points A. with M. = 1,2,3s 
...to only advance a. small step and for each new set of attabhment 
points we will use the results of the previous calculation as the 
zero approximation for the superposition factors: 

New = previous b! and 

1 d'*’ -d* ^ 

hew and previous . dj^ , depending on whether for pro- 

file number 1 the attachment point wis displaced forwards or not. 


At the end of this section we will now discuss how reasonable /24 
zero approximations of the superposition factors can be obtained 
for the first attachment points lying near the trailing edges. 


If the displacement region near the trailing edge has the 
height h which Is simulated by the outflow region (see Figure 4) and 
the average speed of the outflow Is approximately V^, then approxi- 
mately the following relationship has to hold: 

(12a) ^ where 

i 

(12b) - J yield of the 

® 0 i source occupation 

With our selected source occupation, the yield Is 

"Fbr small displacement regions, then from equation (12a) we find 

° (h/t)/o.7. ^j^g estimate h/t = 1.2,...,N, we obtain 

We will select the superposition factor In such a way that the 

resulting overall flow provides for smooth departing flow at the 
trailing edge of each profile 1: 
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“tH 


(aj) (a^) (bj^) 






^ ° ^VtH^ " 

* E *‘*j,i“tHj' ’ = °‘ 

j = 1 

in 


Based on the additio.nal conditions selected for the basic solu- 
tions, all of the are zero, with the exception of 

and . Therefore, equation (l4) gives 

“tH^' = ° 

Using equation (13) we find 


b|»> » -0,6 . ^ 


i = 1,2,... ,N 


/25 


We now have the required approximation values for the first Itera- 
tion. 


2.2.3 Finding separation point 

After calculating the resulting flow with the separation 
regions behind the attachment points *N,mjj 

boundary layer calculations are carried out for the top sides of 
all profiles. For each profile 1 these determine whether the 
boundary layer would separate for this flow ahead of . 

If this Is the case, then for the corresponding profile the point 
6 Is considered as new attachment points for the next cal- 
culation which should lie a small distance ahead of It; otherwise, 
the previous attachment point Is retained. Our attachment points 
are placed along the upper side support points along the profile 
contours In such a manner that when m^^ Is Increased by 1, the attach- 
ment point ’.moves to the next forward support point by a small amount. 

The overall calculation Is repeated for a new set of attachment 
points. This Is repeated until a set Is found 'for which the result- 
ing flow for the first time no longer has any separation ahead of 
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any of the attachment points. This flow Is the result of our cal- 
culation. Since for the attachment points just behind them, separa- 
tion would have resulted ahead of these points now that It no 
longer occurs, one would expect that the last used attachment points 
will agree with the desired and physically correct separation 
point. One has to determine, however, whether one can justify 
always moving forward when finding the correct set of attachment 
points of the separation region (or whether one should stand still). 

We, therefore, have assumed that an Increased separation at one of 
the profiles will also result In an Intensified separation at the 
other profiles (or at least It will stay the same), but will never 
result In a reduced separation. Let us consider the most Important 
case In practice of a slat-main wing — slotted flap combination 
with the following possibilities: /26 

(a) Intensified separation at the flap: 

Because of the reduced flap effect (reduced speed at the rear 
part of the main wing topside) Intensified separation can be expected 
over the main wing. 

(b) Intensified separation at the slat: 

Because of the reduced slat effect (reduced unloading of the 
main wing In the front part, higher suction peak) Intensified separ- 
ation can be expected over the main wing. 

(c) Intensified separation over the main wing: 

The reduced circulation over the main wing again gives a flow 
which has a reduced Inclination but a slower flow In the region of 
the slat and a flow with higher Inclination In the region of the 
(extended) flap. For the flap It Is not certain which effect dom- 
inates but for the flap. Intensified separation can be expected. 


21 



From (b) and (c), we finally find the Indirectly Intensifying 
Influence of a slat separation on the flap separation. 

This means that with some care the assumption made above for 
a tendency of alternating Intensification of the separation regions 
can -be accepted. 

The other question has to be posed: Is the final .solution 

found with the method described the only possible one or only that 
solution with the smallest possible separation regions? 

The question of whether there are physically correct solutions 
with larger separation regions will be discussed In section 3*4 
for several calculation examples. When looking for the separation 
points, we will use the methods discussed above as well as the 
reverse method by starting with very large separation regions. 

The attachment points will be displaced backwards In small steps 
until for the first time' separation occurs ahead of all attachment 
points. In 'this way we will determine the solution with the largest 
possible separation region. We should report here that for all of 
the Investigated cases, the solutions with the smallest possible 
and the largest possible separation regions approximately agree, that 
Is j the solution was unique. 

2.2.4 Closed dead water model for Improving the dead water 
pressure calculation and drag calculation 

In [10] for a single profile It was found that a certain Improve 
ment to the calculated results can be achieved, especially for the 
dead water pressure and for the drag, by transferring from the open 
dead water model to a dead water model with finite length. This Is 
achieved by applying an additional sink at a suitable position behind 
the profile. 



For our profile systems with multiple separations, this pro- 
cedure Improvement Is only partially adapted. In order to avoid 
complication, a refined dead water model Is only used for the pro- 
file of the system which Is the farthest downstream. 

2 . 3 Consideration of separation bubbles 

In addition to the previously treated trailing edge separation 
which expands from the trailing edge with Increasing Inclination 
and has a displacement region which extends above the trailing edge, 
the separation can also occur In the form of so-called separation 
bubbles In which the separation region Is limited by a reattachment 
of a flow along the profile contour. One distinguishes short and 
long separation bubbles. 

Short separation bubbles have a length of a maximum of only a 
few per cent of the profile chord length, are rather thin and have 
no Important displacement effect. They can Influence the flow In 
two ways: First of all, by changing the boundary layer and, there- 

fore, the trailing edge separation. Secondly, by a sudden bursting 
which can generate a large separation region and, therefore, a sub- 
stantial lift loss. Short bubbles are caused by laminar separation 
and turbulent reattachment of the flow. According to experience, 
they occur for medium Reynolds numbers at profiles with moderate 
thickness and slight curvature, along the profile topslc^e near the 
nose. /28 


In [l4] It Is shown for a single profile how the short bubble 
can affect the trailing edge separation approximately. This Is done 
by using empirical criteria for the occurrence and bursting of short 
bubbles, as well as empirical formulas for the bubble length and the 
Initial values of turbulent boundary layer behind the bubbles (detail 
see [l4]). These criterion formulas all contain parameters at the 
laminar separation point ' which are avail- 

able after the conclusion of the laminar boundary layer calculation. 
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The extension of the method could occur alone within the framework 
of the boundary layer calculation: 

During the calculation of the laminar boundary layer, it is 
examined whether transition has occurred. If not, then at the lami- 
nar separation point LA, it is examined whether a short bubble can 
exist. If this is the case, then its length Ig is calculated and 
the turbulent boundary layer calculation is started only at a dis- 
tance Ig behind LA. The momentum and energy loss thickness of LA 
are used, as initial values. Otherwise, the turbulent boundary layer 
calculation is started at the transition point or at LA. If burst- 
ing of the bubble can be expected, then LA is .assumed to be a turbu- 
lent separation point. 

The expanded boundary layer program was used in our pi-oflle 
system calculations. If a short bubble exists, one usv Oly .obtains 
a somewhat reduced trailing edge separation -.and. .;.t-.herefcre.> a eo:;?- 
what larger lift. 

Long bubbles can extend over a large part of the pf'oi'lle 'jhord, 
can be very thick and can have a substantial direct d.lspiacemret 
effect on the external flow. If such bubbles are on the profile 
topside, which especially occurs for thin profiles . or small Reynolds 
numbers, then there Is .a li ft decrease and drag. lncre.ase . . Eaiiable 
calculation methods for such bubbles are not yet; availnble.. For pro- 
file systems especially at higher Reynolds numberti r u-h tepside 
large bubbles probably occur rarely. More often larg, tubules occur 
along the profile bottomslde and when they separate, t.noY are 
triggered at the specified points with large curvature (for example, 
at the slat bottomslde and at the beginning of the flap cavity). 

Such bubbles on the bottomslde can substantially chang; the pressure 
distribution, but usually are not very important for- the overall flow 
and the lift. They can be considered with estimated contour modifl- /29 
cations to sufficient accuracy. We will discuss this in section 3*2. 
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2 . 4 Consideration of compressibility 


High lift systems are being used for low speeds (the Lande- 
Mach numbers of large commercial aircraft are * 0,20)), Because 
of local over-speeds, compressibility Influences can be felt 
already for Incident’ Mach numbers starting with ^ ° • -For a 

pure subsonic flow, because of Increased pressure gradients, the 
trailing edge separation Is somewhat Increased and, therefore, there 
Is a reduction in the lift. Local supersonic regions and compress- 
ion shocks can occur already starting at * o» 2 o 

We will now consider a relatively simple compressibility cor- 
rection and will restrict ourselves to pure subsonic flow. This was 
developed in [ 15 ] and was successfully used for single profiles. 


After obtaining the final solution for incompressible flow with 


separation, the pressure coefficient c^ 
ing to [15] using the approximate formula 


is recalculated accord- 


( 15 ) 


°p ink. = 1 ~ I 1 


2 


Ma^ 


1 . 

2 » p 


to find the pressure coefficients c In the compressible flow. (The 

P 

calculation of c from this equation is done by iteration by solving 

P 

the equation ink."^ ° ° using the Newton method). Using 
the pressure coefficients corrected in this way, the boundary layer 
calculation Is again carried out with consideration of the compress- 
ibility. The search for the physically correct solution Is continued 
if separation now occurs ahead of at least one dead water attachment 
point. 


In order to avoid cases where the given incident Mach number is 
over-critical during each compressibility correction it- Is examined 
whether the pressure coefficient c^ is less than the value /30 

°p ink 
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using equation (15) corresponding to the critical pressure coeffi- 
cient 


( 16 ) 


KMa’ 


I 1 - ( -l-T * -^Ma* ) 

K+1 »C + 1 ® 


ic-1 


anywhere. In this case, the Mach number is reduced until everywhere 

ft 

we have ®p ink. ^ ink. • In the case of a given over-critical 
Mach number, the critical Mach number is determined and the calcu- 
lation is carried out for it. 


3. METHOD TESTING AND IMPROVEMENT TESTS 


In order to test our calculation method, its application limits 
and possible improvements, numerous examples were calculated and the 
theoretical results were compared with measured results as far as 
possible . 


3.1 First results for a profile with Fowler flap 

The first flow calculations were done for the profile GA(W)-1 
with a 30^ deep Fowler flap, for which experimental investigations 
were described in detail in [l6,17]* This Included force measure- 
ments, pressure distribution measurements and flow visualizations 
for various flap positions and angles of attack, but for moderate 
Reynolds numbers (between Re = 2 • 10“ and Re = 3 * 10°). The Mach 
number was = o,i5 . Xn order to fix the laminar-turbulent trans- 
ition, turbulent strips were used for all measurements which were 
attached at 5% profile chord along the topside and the bottomslde 
of the main wing. 

Figure 5 shows the lift curves according to theory and experi- 
ments in the case of deployed flaps with 30^ flap deflection and 
for comparison also for the case of the retracted flap (single pro- 
file). The Reynolds number was 2.2 • 10°. For the calculations, 
it was first assumed that there were no separation bubbles on the 
topside and that the transition occurs at the laminar separation 
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point but at the latest at the location of the turbulent strips, 
and also that the flow was incompressible (Ma^ = 0). For the case /31 
of the deflected flap, for small to moderate angles of attack there 
is a good agreement between the theoretical curve and the experiment- 
al value (symbols). .In particular, it is found that the substantial 
increase in the lift compared with the single profile as well as the 
lift loss (caused by friction and partial separation) as compared 
with ideal flow is well represented by our theory. The maximum lift 
coefficient is substantially Increased by the deployed flap as 
follows ; 

according to measurements from about 1.55 to about 3.6 and 
according to theory from about 1.50 to about 3.8. 

This result is somewhat satisfactory as well. On the other 
hand, in the post stall range, the results are not satisfactory for 
the present profile system. We will discuss this later. 

By consideration of the compressibility (dashed curves), the 
results were somewhat improved but not drastically. 

The following observations are unsatisfactory. 

a) For the single profile at moderate angles of attack, the 
measured lift coefficients are clearly higher than the calculated 
ones. This could be caused by the fact that here short separation 
bubbles occur in practice in spite of the turbulent strips which 
already occurred ahead of the turbulent strips. A second calculation 
where short separation bubbles are allowed resulted in the thin 
curve shown. According to this, such bubbles seem to also exist 

for moderate angles of attack but not in the post-stall range. Our 
empirical criterion for the existence of bubbles also gave short 
bubbles in the post-stall range but is unreliable here. 

b) For the flap profile in the post-stall range, the calculated 
lift coefficients are much higher than the measured ones. In the 
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present example physical phenomena must come Into play, which are 
not completely covered in the theoretical model yet. In order to 
obtain a better Insight Into the physical reality, In Figure 6 we 
show the pressure distribution of the flap profile for a =-10.3® stnd 
12.8^. At a = IO. 3 .J the theory and the measurement have only small /32 
separations over the main wing. At the flap our calculation results 
in moderate trailing edge separation, whereas the measured pressure 
distribution indicates a separation bubble near the nose. Overall, 
the pressure distributions agree for the most part. 

At a = 12.8° on the other hand, the measurement at the flap 
shows a remarkable pressure distribution. It deviates from the 
pressure distribution for the flow without separation (dashed curve) 
as well as from the pressure distribution calculated for trailing 
edge separation (solid curve) to a significant degree. Apparently, 
we have no clear trailing edge separation here but a very strong 
Influence of friction. Perhaps this is a long bubble with adjacent 
very thick boundary layer. Over the main wing the measurement shows 
clear trailing edge separation with substantial dimensions, whereas 
the calculated flow is approximately attached. The unrealistic 
modeling of the flap flow here apparently leads to a false result; 
with respect to the position of the separation point over the main 
wing. 


Clear agreement in the pressure distribution over the main 
wing, on the other hand, occurs when the position of the separation 
point is correctly specified there (solid curve in Figure 7) • For 
the entire flow in this case, the lift coefficient of °a “ 
results which is not very different from the measured one. At the 
flap the calculation gives 50 % trailing edge separation with a 

irj 

For the calculation the computer program was extended so that the 
iteration determination of the separation points could be selected 
for those profiles, for which no fixed specified separation points 
were available. 


28 



pressure distribution which is clearly different from the measured 

t 

one. 


In order to examine to what extent the modeling close to real- 
ity of the flap flow might lead to a better overall solution,' we 

% 

attempted to achieve a pressure distribution over the flap with a 
suitable contour modification which is close to the measured one. 

The result is the dashed curve in Figure 7. Automatically, we found 
a stronger separation and a better pressure distribution over the 
main wing and a lift coefficient of c =2.6 which almost agrees with 
the measured one. This means that we would be in the position of /33 
calculating the overall flow even for substantial separation over 
the main wing (that is, in the post-stall range), if we were only 
able to analyze the flow over the flap better. For this purpose, 
we would have to have more information about how the measured press- 
ure distribution comes about over the flap and we would have to know 
the responsible physical phenomena. The following possibilities 
should be considered: 

a) long bubble over the flap topside 

b) insufficient modeling of the main wing wake and possible contact 
of this wake with the flap boundary layer (wake boundary layer 
confluence) 

c) Influence of separated flow at the flap cavity on the flow around 
the flap. 

Better Information about the true condition can only be obtained 
with very careful experimental work on the flow field in the vicinity 
of the flap. In the present case, we could not progress using 
theory alone. It is possible that the complications are caused for 
the most part by the relatively low Reynolds numbers of the measure- 
ment and the special geometry of our example. For large Reynolds 

7 

nutibers over 10 which are Important for large aircraft, topside 
separation bubbles probably would no longer occur and the friction 
effects would be in general smaller. This fact and somewhat larger 
gap dimensions between the main wing and the flap make a wake- 
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boundary layer confluence less probable. Therefore, we should 
expect that our method in its present form (that is, without con- 
sideration of phenomena mentioned in a) and b)) will give better 
results for large Re and more favorable configurations. 

3.2 Tbe problem of flap cavity flow 

Along tbe bottom side of the main wing at the beginning of the 
flap cavity there is a corner or at least a strong curvature. In 
both cases one would expect separation of the flow. For the re- 
flected flap, the flow which is directed upwards, ahead of the flap 
is probably reattached again. In the flap cavity, therefore, we 
expect a separation bubble. In order to estimate the displacement 
effect of such a bubble, we carried out calculations for two extreme 
assumptions as follows: 

a) for Inflnitestlmally small bubbles (the original small con- 
tour of the flap cavity was used as the profile contour) 

b) for a thick bubble whose contour corresponds to the retracted 
flap (the original profile GA(W)-1 without flap cavity was used as 
the main wing in the calculation) . 

The pressure distributions a = k9,3° for both cases are compared 
in Figure 8 and the measured ones. It can be seen that the bubble 
only has a local influence on the pressure distribution. It can 
also be seen that the simple assumption of a bubble which has the 
shape of the retracted flap already leads to a rather usable correct- 
ion of the pressure distribution along the bottomside of the main 

wing. This assumption was, therefore, retained for all previous 

» 

calculations . 

Regarding the displacement effect, it seems that a very accur- 
ate modeling of the flap cavity bubble does not seem to be necessary. 
_ 

'Using a somewhat shorter bubble, which reattaches at about 5^ 
ahead of the trailing edge and which has a maximum convex curvature 
at x/l=0.75> one obtains a slightly better matching to measured 
results . 


30 



On the other hand, other Influences should not be excluded, for 

9 

example, that the transition zone between the bubble and the exter- 
nal flow extends down to the flap boundary layer and It Is in- 
fluenced (due to increased total pressure and/or the change turbu- 
lence of the incoming flow) . Also there is the possibility of 
unsteady separation in the flap cavity, especially for a continuous 
transition. Extensive experimental work about these questions would 
be very interesting in order to obtain more . knowledge and to Improve 
the theoretical model. 

3 . 3 Testing of various trailing edge conditions for 
calculating the dead water region 

As mentioned in section 2.2, we used the following two condi- 
tions for specifying the circulation and the displacement for each 
profile 

(17 a,b)^ °pA “ °pu °pu “ °po 

A and U were the topside and bottom side separation points and 0 was /35 
the point on the upper boundary of the displacement region. U was 
always assumed at 98.5^ of the profile chord and 0 was assumed at a 
distance ^ from A on the streamline which starts at A. There is a 
certain arbitrariness in the selection of points U and 0. 

Using a few examples, we examined whether by changing the condi- 
tions mentioned above the result of our calculations would be sub- 
stantially Influenced and whether they could then be Improved. We 
were primarily Interested in the correct lift value and also in an 
improved determination of the dead water pressure which was usually 
too high in our previous calculations. The test calculations' were 
first performed for the profile GA(W)-1 (with a retracted flap) . 

.First of all, we varied the position of the lower comparison point U 
as follows (see also Figure 9): 

1 ) u,= Ugy, that is poiht U at 97 % of the profile chord along 
the bottom side profile contour. 
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that is point U at 98.5^ of the profile chord 
along -the bottomside contour 

point U is \% behind the profile along the 
extended profile central line. 

ft 

Figure 9 shows the results of test calculations for a = 8° 

(slight separation) and a = 16° (moderate separation). • First of 
all, one can see that compared with the theory of frictionless flow, 
our theory with friction separation in every case results in a sub- 
stantial improvement. 

The different position of the point U brings about slight 
changes in the results for small separation regions (a = 8°) and 
somewhat larger differences for large separation regions (a = 16°). 
Regarding the pressure distribution outside of the dead water region 
and regarding the lift, best agreement with the measurements occurs 
for ** = '*98,5 . The position of the separation point A is well cal- 
culated, but the value of the dead water pressure is clearly too 
high in all cases. “ = “ioi at a = l6° gives the most favorable dead 
water pressure, but leads to a lift coefficient which is clearly too 
high. In addition, sometimes there are convergence difficulties with 
« = “loi . 

/36 

Additional test calculations of this kind were performed with 
the profile NACA ^l4l5- They confirm the above. The most important 
input data and results are given in the following table. 


2) U £ 

3) U 5 


PROFILE NACA 44l5 Re = 3 * 10'’, Ma^ < 0.15 

(calculation for: incompressible flow, transition in the lam. 

separation point) 

a 

Lift coefficient c according to: 

wlthoilt meS^"'^' theory with friction and separation 

friction from [l8] U e U 97 U = Uba-s U= Uxox 

8° 

1^1° 

1.52 1.20 1.27 1.26 1.25 

2.24 1.40 1.37 1.39 1.45 
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The table again shows a very slight dependence of the result 
on the position of the point U at a = 8® and only a moderate depen- 
dence at a = l4°. The most favorable result is again obtained for 


In this way we calculated the complete lift curves for two 
different Reynolds numbers and these are compared in Figure 10 with 
measurements. 

Even though the lift can be easily calculated (as well as the 
separation point), there is still the problem of the excessively 
high dead water pressure. For large dead water regions, this is 
important for the moment and the drag. For the moment curves shown 
on Figure 10 bottom, we find good agreement between the theory and 
the measurement only up to a a 12 *, , and beyond this the calculated 

moment is too large. 

Our attempts to change the conditions (17a, b) in order to ob- 
tain better results were continued. 

As the next step, we varied the position of the point 0 by 
changing the distance AO. Calculations for the profile NACA ^1^15 
gave the following for Re = 3 * 10^ and a = 1^° . /37 


for 

AO 

• 0 , 9 AH 

1^0 AH 

1,1 Aii 

ineasumiier 

values ' 


> 1,390 

1,392 

1,393 

,1,410 

and. 


w 7* 0,020 

- 0,021 

. -0,021 

-0,070 


It is found that the results are rather insensitive to a moder- 
ate variation in the position ol 0 and that the lift agrees very 
well with measured values for all cases. Therefore, we cannot bring 
about a substantial Improvement in the moment coefficient in this 
way. 
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It Is found that the results are rather Insensitive to a 

» 

moderate variation in the position of 0 and that the lift agrees 
very well with measured values for all cases. Therefore, we cannot 
bring about a substantial improvement in the moment coefficient in 
this way. 

Finally, the conditions (I7a,b) were replaced by 
(18a, b) Cp^ = CpQ and = 0, 

In other words, equal pressure is only assumed for points A and 0 
and the requirement for equal pressure in U is given up and 
replaced by the requirement for smooth departing flow at the trail- 
ing edge. The idea of this is to avoid the coupling of the dead 
water pressure and the relatively high pressure along the profile 
bottomslde. Calculations with these conditions gave the following 
for the profile NACA 4415 at Re = 3 * 10^ and a = 14° 

c„ = 1.440 
a 

and c^ = -0.031. 

Here we find an Improved moment coefficient which is still not 
satisfactory at a somewhat higher lift coefficient. Calculations 
for additional angles of attack showed that the same is true for 
larger angles of attack, but for smaller angles of attack, the con- 
ditions (17a, b) and (l8a,b) give practically the same good results. 

Since we do not have anything better, we will continue to use 
the conditions (17a, b), that is U = U^g ^ and AO = 1.0 • AH. 

The investigations in this section showed us that the lift is • 
usually calculated correctly and it changes only slightly when one 
varies U and 0 moderately. The problem of the excessively high 
dead water pressure and moment coefficient for large angles of 
attack, however, remain. We can probably only find a remedy -by 
using a fine dead water model. The model of DVORAK and MASKEW [183 
seems promising where the boundaries of the displacement region 
are represented by vortex layers. 



3. ^ Various methods for finding separation points and 

t 

Investigations about the uniqueness of the solution 

We will now discuss the question whether the separation points 
which are obtained for the searching method discussed In section 
2 . 2.3 for profile system are the only possible ones. It Is hardly 
possible to give a general proof of this. However, there Is one 
way to clarify the question of a unique solution for a .single Indi- 
vidual case. If one accepts the principle of mutual Intensification 
of separation regions, then the previous method (search method 1) 
for which the search starts with very small separation regions, will 
lead to the solution with the smallest possible separation regions 
compatible with one another. Conversely, one will obtain the solu- 
tion with the largest possible separation regions compatible with 
one another If one starts the search with very large separation 
regions and If displacement points are displaced In very small steps 
backwards until, for the first time, separation is established ahead 
of all attachment points (search method 2). If calculations are 
made according to both search methods and if both result in the 
same solution (approximately), then this solution Is unique. 

Strictly speaking. If a unique solution existed, search method 2 
should lead to a somewhat smaller separation region than search 
method 1 because the attachment points A were displaced In finite 
steps, forwards in the case of search method 1 until they were just 
ahead of the separation points.; For search method 2, they were dis- 
placed backwards until they were all just behind the separation point. 

After an appropriate program extension, calculations were mad^ 
according to both search methods for several examples. Figure 11 
shows the results for our GA(W)-1 profile with 30 % Fowler flap and 
for 30 % flap deflection and a = l6°. We show the positions of the 
attachment point for each Iteration step. As expected for the case 
of a unique solution, the attachment points for; both profiles for 
search method 2 are somewhat behind those found for search method 1. 

A solution with larger dead water regions, therefore, does not /39 

exist . 
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In the same way two further cases with multiple separation 
were Investigated, that 'Is 

the slat-main wing configuration of Figure 13 for 

Re = 3,3 • 10*, Ma^= 0, a = 36‘ (strong Separation over the main 

wing, slight separation over the slat), and 

the three part configuration consisting of slat, main wing arid flap 
of Figure l4 for Re = 1,5 • 10 ’, Ha^= 0 , a = 20 ® (strong flap separ- 

ation, small separations over the slat and the main wing) . 

In both cases we found uniqueness of the solution. 

With the assumption that the solution Is always unique, of 
course only a single search method has to be used. Search method 1 
should then be preferred because, compared with search method 2, 

It usually requires less computation time. Also, search method 2 
sometimes has convergence difficulties In the Iteration calculation 
of the superposition factors. 

In the calculations with slats It was found that the principle 

of mutual Intensification of separation regions Is only partially 

applicable for the slat. The slat separation does Influence the 

separation over the main wing but conversely, for substantial main 

* 

wing separation there is a somewhat reduced slat separation . From 
this the requirement arises for extending the search method in such 
a manner that the dead water points can migrate not only forwards 
during the search process, but If necessary, also backwards. This 
was brought about by using search method 2 after using search method 
1 with the found attachment points. Starting with Its result, 
finally search method 1 is used again. In this way by searching 
forward?, backwards and forwards, the correction to those separation 
regions 1?^ brought about which would be too large If only search 
method 1 or 2 were used individually. 

*]^xplanatlon : Reduced circulation over main wing results In a 

smaller effective angle of attack of slat. 
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For our three part configuration. Figure 12 shows the itera- 
tion for the combined search method. 


The following table gives the positions of the dead water 
points along the various wing parts according to the various search 
methods . 



slat 

main wing 

flap 



X 

II 

-0.025 

0.862 

1.037 

according to 
method 1 

search 


-0.020 

0.872 

1.037 

according to 
search method 

combined 

(1-2-1) 


-0.000 

0.880 

1.052 

according to 
method 2 

search 


As expected, search method 2 always gives values which are 
somewhat too large (that is, separation regions are too small). 
There are differences between the results of search method 1 and 
the combined search method only for the slat and the main wing and 
they are not large. 

For safety, in the following we will always use the combined 
search method. There is a possibility of reducing the calculation 
effort as will now be discussed. 

In Figure 12, it becomes apparent that the steps should be 
substantially increased during phase 1 of the searching process. 
This probably is acceptable for the combined search method as well 
because in this way the true separation point is clearly exceeded 
and because this can again be connected by the backward searching 
in phase 2. On the other hand, it has to be considered that for 
large steps, the superposition factors b and d will change drasti- 
cally from step to step so that for the inner iteration (determina- 
tion of the superposition factors) the zero approximation values 
will be not as good and, therefore, one would expect higher compu- 
tational times and even convergence difficulties. It does not 
seem possible to give general Information about the optimum step 
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size with respect to the overall computation time, but by carefully 
adjusting the step size. during search phase 1, a’reductlon In the 
overall computational time can probably be brought about. There- 
fore, for the present example, we carried out sample calculations 
with the double and triple step size. Both calculations gave the 
same results as for simple step size (without convergence diffi- 
culties). The overall computational time (using the IBM 370/158 
computer) were the following: 
for the combined search method 


with simple step size 

19 

minutes 

CPU time 

with double step size 

17.5 

It 

II 

with triple step size 

16.6 

II 

II 

for search method 1 

18 

II 

II 

for search method 2 

26 

II 

II 


/4l 


Accordingly, the combined search time with the triple step size 
is the most favorable. A further increase in the step size does not 
bring about any substantial calculation advantages but increases the 
danger of convergence difficulties for the Inner iteration. The 
possible calculation time gained by further Improvement of the 
search method is too limited because most of the total calculation 
time (about l4 minutes in the present case) Is required outside of 
the search process, especially for the calculation of the basic solu- 
tions. 


For the further calculation we use the combined search method 
with the triple step size during search phase 1 (and simple step 
size during the other phases). 

At this point we would again like to point out that in the 
case Ma^g / 0 the compressibility correction (according to section 
2.4) is used in the last two search phases whereas the (rough) 
search phase I'l.s only carried out for incompressible flow in order 
to save time. 
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3 . 5 Calculation results for profile with slat 


For the two part configuration consisting of slat and main 
wing Investigated experimentally In [20], Figure 13 shows the lift 
curve according to theory measurements for Re = 3*4 * 10^ and 
Ma^ = 0.l6 for the two gap widths s. For comparison, we also show 
the lift curve for the fast flight profile NACA 64-210 without 
slat. The measured and calculated maximum lift coefficients 
agree In all cases very well. Both the substantial lift gain with /4 2 
the slat as well as the slight displacement of the lift maximum with 
the change In the gap are well predicted. In the post-stall range, 
on the other band, the calculated lift coefficients are always some- 
what higher than the measured ones. 

The lift maximum Is exceeded also when a slat exists due to 
advancing separation over the main wing. The small separation over 
the slat only starts to expand above a s 35“ , Por the calculations. 

It was, therefore, assumed that no bubbles exist and that the trans- 
ition occurs at the laminar separation point. For the main wing, 
we used = 120 target points and for the slat we used M 2 = 60 for 
target points. The calculation time on the IBM 370/158 was 7-11 
minutes per angle of attack, depending on angle of attack (and, there- 
fore, separation regions). 

3 . 6 Calculation results for a three-part profile system 

For a three-part high lift configuration consisting of slat, 
main wing and Fowler flap designed at the DFVLR Institute for Design 
Aerodynamics, pressure distributions and lift curves were calculated. 
Measurements for these configurations are In preparation but experi- 
mental results are not yet available. 

Calculations were first done for rather large Reynolds numbers 
(Re = 1.5 * 10 ), as might occur during takeoff and landing of large 
cruise aircraft. Figure l4 shows the configuration and the calcu- 
lated pressure distribution for a = 0° and a = 20° (corresponds 
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approximately to c ). For a = 0°, there was a substantial sep- 

aration on the topside only over the flap. But the following bound- 
ary layer calculation which is also carried out for the bottomslde 
of the profile, shows turbulent separations at the bottomsldes of 
the slat and the main wing. With the assumption that in this way 
long bubbles are created in the slat bottomslde and in the slat 
cavity, the calculation was again repeated with the corresponding 
modified profile contours. The resulting pressure distribution is 
also shown in Figure l4 and shows the displacement effect of such 
bubbles on the bottomslde. Substantial changes in the pressure dls- /^3 
tribution only occurred near the bubbles. Over the flap there was 

only a slight enlargement of the suction rpeak. The separation points 

* 

and the total lift coefficient practically did not change . 

For a = 20° on the topside, there was a strong separation at the 
flap and In addition, no separations over the main wing and the flap, 
whereas on the bottomslde there was no separation. As Figure 15 (top) 
shows, as the angle of attack is Increased further (a > 20°), there 
is an intensification of flap separation and a rapid Increase in the 
separation over the main wing, and therefore, a reduction in the 
lift coefficient. The separation over the slat first Increases 
slightly and then decreases slightly again for larger angles of 
attack. Figure 15 shows the various lift curves in the lower part. 

7 

Except for Re = 1.5 * 10 , the lift curve was also calculated for 
Reynolds numbers which can be achieved In.i. the wind tunnel (Re = 

6 • 10^). There is a strong influence of the Reynolds number on the 
maximum lift of these high lift configurations. Accordingly, for 
the full scale version, one would expect about 13^ higher lift coeff- 
icient than for model measurements in the wind tunnel, provided that 
in both cases there were no physical phenomena which were not 

*■5 

This does not mean that such bubbles could not change the overall 
flow in another way than by their displacement effect. In parti- 
cular, for small gap width it is possible that the transition 
zone between the bubble and the external flow would have an effect 
on the profile behind it because of the reduced energy and the 
changed turbulence in the boundary layer and in this way, the separ- 
ation in the overall flow would be changed. 
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included in the calculation, such as bubble effects and wake-bound- 
ary layer confluence. 

Finally, the calculation was made with consideration of com- 
pressibility (dashed. curve) that is for Ma^ = 0.19. This is approx- 
imately the (lower) critical Mach number at the lift maximum for 

7 

Re = 1.5 • 10 . Because of the influence of compressibility, the 
lift maximum is displaced towards smaller angles of attack and is 
substantially reduced. 

For these calculations, it was always assumed that the laminar 
turbulence transition occurs at the laminar separation point. 60 
target points were used for the flap and 100 target points were used 
for the main wing. The computation time on the IBM 370/158 was 15- 
20 minutes for a set of parameters (a. Re, Ma^) . /44 

i|. SUMMARY AND OUTLOOK 


In numerical calculation method [6] for steady Incompressible 
flow over profile systems with Simple trailing edge , separation was 
extended to cases where there are several separation regions over, 
the various parts of the profile system. It is only by describing 
this multiple separation that it becomes possible to calculate the 
lift curve in the vicinity of the maximum lift and in the post-stall 
range for mechanical high lift systems. 

In the extended calculation method, in addition to multiple 
trailing edge separation, the influences of bubbles and compressi- 
bility were Included for subsonic flows, approximately. The expanded 
method was programmed in FORTRAN and was tested with various examples. 
In all investigations, plausible unique solutions resulted. Compari- 
son between theoretical results in measurements showed for the most 
part satisfactory agreement. 
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It was also found that the calculation possibilities are sub- 
ject to some important restrictions: 

a) steady plane flow is assumed. Unsteady effects for wing 
oscillations, gusts, acceleration, pulsating separation regidns, 
etc., as well as thr^ee-dimensional effects over wings with finite 
stand are not covered. 

b) It is assumed that the wakes of upstream wing parts (wake- 
boundary layer confluence) do not substantially influence the bound- 
ary layer directly. In practice, such influences often exist for 
moderate Reynolds numbers and small gap width. 

c) Compressibility influences are only considered for pure sub- 

sonic flow and- only approximately. In practice, local supersonic 
regions could certainly occjir for high lift configurations. /4 5 

d) Separation bubbles are only considered incompletely. On the 
profile topside, only short bubbles are allowed and these are only 
allowed to a limited extent and only as regarding their effects on the 
trailing edge separation. The empirical criteria present for the 
occurrence and bursting of such bubbles are incomplete and not very 
reliable . 

Bubbles along the profile bottomslde can only be roughly covered 
by estimated contour modifications with respect to their displacement 
effect. In practice, separation bubbles of all kinds occur, espe- 
cially for not very large Reynolds numbers, and can substantially 
Influence the overall flow under some conditions. 

e) Our potential theory dead water model satisfies the condi- 
tion of constant pressure only approximately (at three points) and 
according to experience, results in dead water pressure which is 
somewhat too high. For larger separation regions in this way, the 
moment and the drag are inaccurately calculated. 

In addition,- for separation over the main wing of a slotted 
flap configuration (section 3.1), a modeling which is too rough for 
the main wing wake could be the reason for unsatisfactory descrip- 
tion of the flap flow. 

f) The computer program can only be applied for trained persons 
in its present form because of continuous Improvements. It is not 



yet optimized and it is insufficiently documented. It still does 
not contain the ground effect. 

Prom all this, we have the following stimuli for future work: 

- Improvement, completion and documentation of the computer 
program 

- development of a refined dead water model for strong and 
multiple separation 

- better coverage of all bubble effects 

- better coverage of compressibility Influences, if possible 
for transonic flows 

- consideration of wake, boundary layer confluence 

- development of methods for three-dimensional and unsteady 
flows with separation 

This further development cannot occur purely theoretically and 
will require substantial support of experimental work. 
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Figure 6 
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Figure 8 
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Figure 9 . 
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Figure 10. 









Figure 11 


GA(W)-1 with 30% Fcwler flap;iF = 30*’ 
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